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Introduction 



A soliton is a localized lump (or string or wall etc) of energy, which can move without 
distortion, dispersion or dissipation, and which is stable under perturbations (and 
collisions with other solitons). The word was coined by Zabusky and Kruskal in 
1965, to describe a solitary wave with particle-like properties (as in electron, proton 
etc). Solitons are relevant to numerous areas of physics — condensed- matter, cos- 
mology, fluids /plasmas, biophysics (eg DNA), nuclear physics, high-energy physics 
etc. Mathematically, they are modelled as solutions of appropriate partial differential 
equations. 

Systems which admit solitons may be classified according to the mechanism by 
which stability is ensured. Such mechanisms include complete integrability, nontrivial 
topology plus dynamical balancing, and Q-balls/breathers. 

Sometimes the term 'soliton' is used in a restricted sense, to refer to stable lo- 
calized lumps which have purely elastic interactions: solitons which collide without 
any radiation being emitted. This is possible only in very special systems, namely 
those that are completely-integrable. For these systems, soliton stability (and the 
elasticity of collisions) arise from a number of characteristic properties, including: a 
precise balance between dispersion and nonlinearity, solvability by the inverse scat- 
tering transform from linear data, infinitely many conserved quantities, a Lax for- 
mulation (associated linear problem), and Backhand transformations. Examples of 
such integrable soliton systems are the sine-Gordon, Korteweg-deVries and nonlinear 
Schrodinger equations. 

The category of topological solitons is the most varied, and includes such examples 
as kinks, vortices, monopoles, Skyrmions and instantons. The requirement for these 
of 'dynamical balancing' can be understood in terms of Derrick's theorem, which 
provides necessary conditions for a classical field theory to admit static localized 
solutions. The Derrick argument involves studying what happens to the energy of 
a field when one changes the scale of space. If one has a scalar field (or multiplet 
of scalar fields) 4>, and/or a gauge field F^, then the static energy E is the sum of 
terms such as 

Eo = J V{<j>) d n x, E d = J T d {D j( t>) d n x, E F = J F jk F jk d n x, 

where each integral is over (ra-dimensional) space R™, Dj(f> denotes the covariant spa- 
tial derivative of <p, and T d (£j) is a real- valued polynomial of degree d. In particular, 
for example, we could have T2{Dj<p) = (Dj(f))(Dj(p): the standard gradient term. 
Under the dilation x J ^ \x 3 , these functionals transform as 

E ' ^ \~ n E , E d ' ^ X d ~ n E d , E F ' ^ X 4 - n E F . 
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In order to have a static solution (critical point of the static energy functional), 
one needs to have a zero exponent on A, and/or a balance between positive and 
negative exponents. A negative exponent indicates a compressing force (tending to 
implode a localized lump), whereas a positive exponent indicates an expanding force; 
so to have a static lump solution, these two forces have to balance each other. For 
n = 1, a system involving only a scalar field, with terms of the form Eq and E2, 
can admit static solitons — for example kinks; the scaling argument implies a Virial 
Theorem, which in this case says that Eq = E2. For n = 2, one can have a scalar 
system with only E2, since in this case the relevant exponent is zero — for example 
the 2-dimensional sigma-model. Another n = 2 example is that of vortices in the 
Abelian-Higgs model, where the energy contains terms Eq, E2 and Ep. For n = 3, 
interesting systems have E2 together with either E4 (for example Skyrmions) or Ep 
(for example monopoles). An Eq term is optional in these cases; its presence affects, 
in particular, the long-range properties of the solitons. For N = 4, one can have 
instantons in a pure gauge theory (term Ep only). 

It should be noted that if there are no restrictions on the fields (ft and Aj (such 
as arise, for example, from non-trivial topology), then there is a more obvious mode 
of instability, which will inevitably be present: (j) 1— > /j,(f> and/or Aj ^ fJ-Aj, where 
< n < 1. In other words, the fields can simply be scaled away altogether, so that 
the height of the soliton (and its energy) go smoothly to zero. This can be prevented 
by non-trivial topology. 

Another way of preventing solitons from shrinking is to allow the field to have 
some 'internal' time-dependence, so that it is stationary rather than static. For 
example, one could allow the complex scalar field 4> to have the form 4> = V>exp(iwt), 
where ip is independent of time t. This leads to something like a centrifugal force, 
which can have a stabilizing effect in the absence of Skyrme or magnetic terms. The 
corresponding solitons are Q-balls. 



Kinks and Breathers 

The simplest topological solitons are kinks, in systems involving a real-valued scalar 
field (j)(x) in one spatial dimension. The dynamics is governed by the Lagrangian 
density 



(for - (fo) 2 - w(<t>y 



where W(fo) is some (fixed) smooth function. The system can admit kinks if W(fo) 
has at least two zeros, for example W(A) = W(B) = with W(fo) > for A < <j> < B. 
Two well-known systems are: sine-Gordon, where W(fo) = 2 sin(</>/2), A = 0, B = 2ir; 
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and phi-four, where W(4>) = 1 — cf) 2 , A = — 1, B = 1. The corresponding field 
equations are the Euler-Lagrange equations for £; for example, the sine-Gordon 
equation is 

<Ptt ~ 4>xx + smcp = 0. (1) 

Configurations satisfying the boundary conditions 4> ^ A as x ^ — oo and 4> — ► B 
as x — ► oo are called kinks (and the corresponding ones with x = oo and x = — oo 
interchanged are antikinks). For kink (or antikink) configurations, there is a lower 
bound, called the Bogomol'nyi bound, on the static energy E[<j>\\ for kink boundary 
conditions, we have 

i roo 

E[<f>] = -z / U x f + W{ct>f]dx 

£ J— oo 

/oo 
W(<p)(t) x dx 
-oo 

> f B w{4>)d^ 

J A 



oo 
oo 



— oo 



with equality if and only if the Bogomol'nyi equation 

| = W W (2, 

is satisfied. A static solution of the Bogomol'nyi equation is a kink solution — it is 
a static minimum of the energy functional in the kink sector. For example, for the 
sine-Gordon system we get E[<f>] > 8, with equality for the sine-Gordon kink 

4>(x) = 4 tan -1 exp(x — xq); 

while for the phi-four system we get E[<j)\ > 4/3, with equality for the phi-four kink 

<t>{x) = tanh(x — Xo). 

These kinks are stable topological solitons; the non-trivial topology corresponds 
to the fact that the boundary value of cj)(t, x) at x = oo is different from that at 
x = — oo. With trivial boundary conditions (say ^Aasn ±oo) stable static 
solitons are unlikely to exist, but solitons with periodic time-dependence (which is 
this context are called breathers) may exist. For example, the sine-Gordon equation 
and the non-linear Schrodinger equation both admit breathers — but these owe their 
existence to complete integrability. By contrast, the phi-four system (which is not 
integrable) does not admit breathers; a collision between a phi-four kink and an 
antikink (with suitable impact speed) produces a long-lived state which looks like a 
breather, but eventually decays into radiation. 
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In lattice systems, however, breathers are more generic. In a 1-dimensional lattice 
system, the continuous space R is replaced by the lattice Z, so <j)(t,x) is replaced by 
(j) n {t) where n G Z. The Lagrangian is 

L = \ E [^nf ~ h~ 2 (<f>n+l ~ Kf ~ W(<P n )] ; 

n 

here h is a positive parameter, corresponding to the dimensionless ratio between the 
lattice spacing and the size of a kink. The continuum limit is h — ► 0. This system 
admits kink solutions as in the continuum case; and for h large enough, it admits 
breathers as well, but these disappear as h becomes small. 

Interpreted in three dimensions, the kink becomes a domain wall separating two 
regions in which the order parameter <f> takes on distinct values; this has applications 
in such diverse areas as cosmology and condensed-matter physics. 

Sigma- Mo dels and Skyrmions 

In a sigma-model or Skyrme system, the field is a map <fi from space-time to a 
Riemannian manifold M; usually M is taken to be a Lie group or a symmetric 
space. The energy density of a static field can be constructed as follows (the Lorentz- 
invariant extension of this gives a relativistic Lagrangian for fields on space-time). 
Let <f) a be local coordinates on the m-dimensional manifold M, let h a b denote the 
metric of M, and let x 3 denote the spatial coordinates on space R n . Define an m x m 
matrix D by 

D a b = (dj^haddjtf), 

where dj denotes derivatives with respect to the x 3 . Then the invariants £2 = tr(D) = 
\dj(f> a \ 2 and £4 = \ [(tr D) 2 — tr(D 2 )] can be terms in the energy density, as well as 
a zeroth-order term £q = V((p a ) not involving derivatives of (ft. A term of the form 
£4 is called a Skyrme term. 

The boundary condition on field configurations is that (f> tends to some constant 
value 4>o € M as \x\ — > 00 in R n . From the topological point of view, this compactifies 
R n to S n . In other words, (j) extends to a map from S n to M; and such maps are 
classified topologically by the homotopy group 7r ra (M). For topological solitons to 
exist, this group has to be non-trivial. 

In one spatial dimension (n = 1) with M = S 1 (say), the expression £4 is iden- 
tically zero, and we just have kink-type systems such as sine-Gordon. The simplest 
two-dimensional example (n = 2) is the 0(3) sigma model, which has M = S 2 with 
its standard metric. In this system, the field is often expressed as a unit 3-vector field 



5 



4> = (0 1 ,0 2 ,0 3 ), with £ 2 = (dj(p) ■ {dj4>). Here the configurations are classified topo- 
logically by their degree (or winding number, or topological charge) N G ir 2 (S 2 ) = Z, 
which equals 

N = — [ $-di$xd 2 $dx 1 dx 2 . 
Air J 

It is often convenient, instead of <p, to use a single complex- valued function W related 
to 4> by the stereographic projection W = ((f) 1 + i(f> 2 )/(l — (f> 3 ). In terms of W, the 
formula for the degree N is 



j_ r w lW2 -w 2Wl x 2 

2vri (l + II^I 2 ) 2 ' 



(1 + |H^| 2 ) 2 
and the static energy is (with z = x l +\x 2 ) 

E - I £ 2 d 2 x 

\W z \ 2 + \W-\ 2 , 2 
■ d x 



(l + \W\ 2 ) 2 



(1 + |H'| 2 )' 2 J (1 + IW'l 2 ) 2 

From this, one sees that E satisfies the Bogomol'nyi bound E > 8ttN; and that 
minimal-energy solutions correspond to solutions of the Cauchy-Riemann equations 
= 0. To have finite energy, W(z) has to be a rational function, and so solutions 
with with winding number N correspond to rational meromorphic functions W(z), 
of degree |7V|. (If TV < 0, then W is a rational function of z.) The energy is scale- 
invariant (conformally-invariant), and consequently these solutions are not solitons 
- they are not quite stable, since their size is not fixed. Adding terms £4 and £q to 
the energy density fixes the soliton size, and the resulting two-dimensional Skyrme 
systems admit true topological solitons. 

The three-dimensional case (n = 3), with M being a simple Lie group, is the 
original Skyrme model of nuclear physics. If M = SU(2), then the integer iV G 
7T3(SU(2)) = Z is interpreted as the baryon number. The (quantum) excitations of 
the (/>-field correspond to the pions, whereas the (semi-classical) solitons correspond 
to the nucleons. This model emerges as an effective theory of QCD, in the limit 
where the number of colours is large. If we express the field as a function U(x 3 ) 
taking values in a Lie group, then Lj = U~ l djU takes values in the corresponding 
Lie algebra, and £ 2 and £4 take the form 

1 
2 



£ 2 = --tr(LjLj) 



£4 = -^tr {[Lj,L k ][Lj,L k ]). 
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The static energy density in the basic Skyrme system is the sum of these two terms. 
The static energy satisfies a Bogomol'nyi bound E > 127r 2 |iV|, and it is believed 
that stable solitons (Skyrmions) exist for each value of N . Classical Skyrmions have 
been investigated numerically; and, for values of N up to about 25, they turn out 
to resemble polyhedral shells. Comparison with nucleon phenomenology requires 
semi-classical quantization, and this leads to results which are at least qualitatively 
correct. 

A variant of the Skyrme model is the Skyrme-Faddeev system, which has n = 3 
and M = S 2 ] the solitons in this case resemble loops which can be linked or knotted, 
and which are classified by their Hopf number N S ^(S 2 ). In this case, the energy 
satisfies a lower bound of the form E > ciV 3 / 4 . Numerical experiments indicate that 
for each N, there is a minimal-energy solution with Hopf number N, and with energy 
close to this topological lower bound. 

Abelian-Higgs Vortices. 

Vortices live in two spatial dimensions, and viewed in three dimensions are string-like; 
two applications are as cosmic strings and as magnetic flux tubes in superconductors. 
They occur as static topological solitons in the the abelian Higgs model (or Ginzberg- 
Landau model), and involve a magnetic field B = d\A2 — 82 A%, coupled to a complex 
scalar field (j>, on the plane R 2 . The energy density is 

£ = \{D^){DJ4>) + \b 2 + ~A(1 - H 2 ) 2 , (3) 

where Dj(j) := dj<j) — iAjcp, and where A is a positive constant. The boundary condi- 
tions are 

D^ = 0, B = 0, |0| = 1 (4) 

as r —* 00. If we think of a very large circle C on R 2 , so that @ holds on C, then 
(j)\ c is a map from the circle C to the circle of unit radius in the complex plane, and 
therefore it has an integer winding number N. Thus configurations are labelled by 
this vortex number N. 

Note that if £ vanishes, then B = and \4>\ = 1: the gauge symmetry is spon- 
taneously broken, and the photon 'acquires a mass': this is a standard example of 
spontaneous symmetry breaking. 

The total magnetic flux J Bd 2 x equals 2-7riV; a proof of this is as follows. Let 
9 be the usual polar coordinate around C. Because \(f)\ = 1 on C, we can write 
<p = exp[i/(0)] for some function /; this / need not be single- valued, but must satisfy 
,f(27r) — /(0) = 2ttN with N being an integer (in order that <f> be single- valued) . In 
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fact, this defines the winding number. Now since Dj(fr = dj<p — iAj<p = on C, we 
have 



If A = 1, then the total energy E = J £d 2 x satisfies the Bogomol'nyi bound 
E > irN; and E = irN if and only if a set of partial differential equations (the 
Bogomol'nyi equations) are satisfied. Since like charges repel, the magnetic force 
between vortices is repulsive. But there is also a force from the Higgs field, and 
this is attractive. The balance between the two forces is determined by A: if A > 1, 
then vortices repel each other; while if A < 1, then vortices attract. In the critical 
case A = 1, the force between vortices is exactly balanced, and there exist static 
multi-vortex solutions. In fact, one has the following: given N points in the plane, 
there exists an A^-vortex solution of the Bogomol'nyi equations (and hence of the 
full field equations) with 4> vanishing at the chosen points (and nowhere else). All 
static solutions are of this form. These solutions cannot, however, be written down 
explicitly in terms of elementary functions (except of course for N = 0) . 



The abelian Higgs model does not admit 3-dimensional solitons, but a non-abelian 
generalization does — such non-abelian Higgs solitons are called magnetic monopoles. 
The field content, in the simplest version, is as follows. First, there is a gauge (Yang- 
Mills) field F^ v , with gauge potential A^, and with the gauge group being a simple 
Lie group G. Secondly, there is a Higgs scalar field cf>, transforming under the adjoint 
representation of G (thus <f> takes values in the Lie algebra of G). For simplicity, take 
G to be SU(2) in what follows. So we may write = iA®a a , F^ = iF^ v a a and 
4> = i(f) a <Ja, where a a are the Pauli matrices. The energy of static (<9o</> = = doAj), 
purely-magnetic (Aq = 0) configurations is 



where Bj = is the magnetic field. The boundary conditions are Bj^O and 

cj) a (j) a —> 1 as r — > oo; so <j) restricted to a large spatial 2-sphere becomes a map from 



Aj = —i<fi 1 dj(f> = djf 



on C. So, using Stokes' theorem, we get 




2-kN . 



Monopoles 
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S 2 to the unit 2-sphere in the Lie algebra su(2), and as such it has a degree N £ Z. 
An analytic expression for N is 

B^(Dj(j)) a Sx = 2itN. (5) 

At long range, the field resembles an isolated magnetic pole (a Dirac magnetic 
monopole), with magnetic charge 2tvN. Asymptotically, the SU(2) gauge symme- 
try is spontaneously broken to U(l), which is interpreted as the electromagnetic 
gauge group. 

In 1974, it was observed pointed out that this system admits a smooth, finite- 
energy, stable, spherically-symmetric N = 1 solution — this is the 't Hooft-Polyakov 
monopole. There is a Bogomol'nyi lower bound on the energy E: from < (B + 
Dtp) 2 = B 2 + (Dcp) 2 + 2B ■ D(p we get 

E > 2ttN + j -A(l - <f<f) 2 d 3 x, (6) 

where use has been made of (jSJ . The inequality © is saturated if and only if we go 
to the Prasad-Sommerfield limit A = 0, and the Bogomol'nyi equations 

{D j( j>f = ~B* (7) 

hold. The corresponding solitons are called BPS monopoles. 

The Bogomol'nyi equations (J2J), together with the boundary conditions described 
above, form a completely-integrable elliptic system of partial differential equations. 
For any positive integer N, the space of BPS monopoles of charge N, with gauge 
freedom factored out, is parametrized by a (4N — l)-dimensional manifold A4n- 
This is the moduli space of N monopoles. Roughly speaking, each monopole has 
a position in space (3 parameters) plus a phase (1 parameter), making a total of 
4|iV| parameters; an overall phase can be removed by a gauge transformation, leaving 
4| N\ — 1 parameters. In fact, it is often useful to retain the overall phase, and to work 
with the corresponding 4 \N\ -dimensional manifold Mn- This manifold has a natural 
metric, which corresponds to the expression for the kinetic energy of the system. A 
point in A4 at represents an iV-monopole configuration, and the slow-motion dynamics 
of N monopoles corresponds to geodesies on A4n', this is the geodesic approximation 
of monopole dynamics. 

The N = 1 monopole is spherically-symmetric, and the corresponding fields take 
a simple form; for example, the Higgs field of a 1-monopole located at r = is 

"coth(2r) 1 " 
r 2r 2 



x 



For N > 1, the expressions tend to be less explicit; but monopole solutions can 
nevertheless be characterized in a fairly complete way. The Bogomol'nyi equations 
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(J7J) are a dimensional reduction of the self-dual Yang- Mills equations in R 4 , and BPS 
monopoles correspond to holomorphic vector bundles over a certain two-dimensional 
complex manifold ('mini-twistor space'). This leads to various other characterizations 
of monopole solutions, for example in terms of certain curves ('spectral curves') on 
mini-twistor space, and in terms of solutions of a set of ordinary differential equations 
called the Nahm equations. Having all these descriptions enables one to deduce much 
about the monopole moduli space, and to characterize many monopole solutions. 
In particular, there are explicit solutions of the Nahm equations involving elliptic 
functions, which correspond to monopoles with certain discrete symmetries, such as 
a 3-monopole with tetrahedral symmetry, and a 4-monopole with the appearance 
and symmetries of a cube. 



Yang-Mills Instantons 

Consider gauge fields in 4-dimensional Euclidean space R 4 , with gauge group G. For 
simplicity in what follows, take G to be SU(2); one can extend much of the structure 
to more general groups, for example the simple Lie groups. Let and F^ v denote 
the gauge potential and gauge field. The Yang-Mills action is 

S = ~jtr{F lxv F lxu )d A x, (8) 

where we assume a boundary condition, at infinity in R 4 , such that this integral con- 
verges. The Euler-Lagrange equations which describe critical points of the functional 
S are the Yang-Mills equations 

DpF^ = 0. (9) 

Finite-action Ynag- Mills fields are called instantons. The Euclidean action (JSJ) is used 
in the path-integral approach to quantum gauge field theory, and therefore instantons 
are crucial in understanding the path integral. 
The dual of the field tensor F^ u is 

the gauge field is self-dual if *F fMU = F^ u , and anti-self-dual if *F flu = —F^ v . In view 
of the Bianchi identity * F^ u = 0, any self-dual or anti-self-dual gauge field is 
automatically a solution of the Yang-Mills equations @ • This fact also follows from 
the discussion below, where we see that self-dual instantons give local minima of the 
action. 

The Yang-Mills action (and Yang- Mills equations) are conformally invariant; and 
any finite-action solution of the Yang- Mills equations on R 4 extends smoothly to the 
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conformal compactification S 4 . Gauge fields on S 4 , with gauge group SU(2), are 
classified topologically by an integer N, namely the second Chern number 

N = C 2 = -g^2 / tr ( F »» * dAx - ( 10 ) 

From (jHJ) and (JTU|) we get a topological lower bound on the action, as follows: 

< - J tr (Ffu, - *F ia ,) (F„ v - *F„ 
= 8S-Wtt 2 N: 



/a/ j d x 



and so S > 2ir 2 N, with equality if and only if the field is self-dual. If N < 0, we 
get S > 27r 2 |iV|, with equality if and only if F is anti-self-dual. So the self-dual (or 
anti-self-dual) fields minimize the action in each topological class. 

For the remainder of this section, we restrict to self-dual instantons with instanton 
number N > 0. The space (moduli space) of such instantons, with gauge-equivalence 
factored out, is a (8N — 3)-dimensional real manifold. In principle, all these gauge 
fields can be constructed using algebraic- geometry (twistor) methods: instantons 
correspond to holomorphic vector bundles over complex projective 3-space (twistor 
space). One large class of solutions which can be written out explicitly is as follows; 
for N = 1 and N = 2 it gives all instantons, while for N > 3 it gives a (5N + 4)- 
dimensional subfamily of the full (8N — 3)-dimensional solution space. The gauge 
potentials in this class have the form 

= ia^dv log <j), (11) 

where the a^u are constant matrices (antisymmetric in fiv) defined in terms of the 
Pauli matrices a a by 



010 


= 023 = 


2 ai 






1 


020 


= 031 = 


— (Jo 

2 






1 


030 


= 012 = 


-03 
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The real- valued function <j> = <j)(x^) is a solution of the 4-dimensional Laplace equa- 
tion given by 



N X 



(„,„ _ ,./< H ,-// _ ,-/' 1 



where the are N+l distinct points in R 4 , and the are iV+1 positive constants: a 
total of 5^+5 parameters. It is clear from (jll|) that the overall scale of cj) is irrelevant, 
leaving a (5N + 4)-parameter family. For = 1 and N = 2, symmetries reduce the 
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parameter count further, to 5 and 13 respectively. Although (p has poles at the points 
x = Xk, the gauge potentials are smooth (possibly after a gauge transformation). 

Finally, it is worth noting that (as one might expect) there is a gravitational 
analogue of the gauge-theoretic structures described here. In other words, one has 
self-dual gravitational instantons — these are 4-dimensional Riemannian spaces for 
which the conformal-curvature tensor (the Weyl tensor) is self-dual, and the Ricci 
tensor satisfies Einstein's equations R^ u = Ag^. As before, such spaces can be 
constructed using a twistor-geometrical correspondence. 

Q-Balls 

A Q-ball (or nontopological soliton) is a soliton which has a periodic time-dependence 
in a degree of freedom which corresponds to a global symmetry. The simplest class 
of Q-ball systems involves a complex scalar field 0, with an invariance under the 
constant phase transformation cf> i— ► e t8 4>; the Q-balls are soliton solutions of the form 

0(t,x) = e^V(x), (12) 

where ^>(x) is a complex scalar field depending only on the spatial variables x. The 
best-known case is the 1-soliton solution 

(j>(t,x) = aV2 exp(ia 2 t) sech(ax) 

of the nonlinear Schrodinger equation ifa + <fi xx + (f>\4>\ 2 = 0. 

More generally, consider a system (in n spatial dimensions) with Lagrangian 

C=\{d^){d^)-U[\ci>\), 

where (f){x^) is a complex- valued field. Associated with the global phase symmetry 
is the conserved Noether charge Q = J lm(</>0t) d n x. Minimizing the energy of a 
configuration subject to Q being fixed implies that (j) has the form (|12|). Without 
loss of generality, we may take u > 0. Note that Q = ul, where I = J |^| 2 d n x. The 
energy of a configuration of the form ()12|) is E = E q + + E p , where 

E l = \J |^| 2 cPx, 
Ek = \lu 2 = ^Q 2 /I, 
E p = Ju(\^\)d n x. 

Let us take U(0) = = U'(0), with the field satisfying the boundary condition ip — ► 
as r — > oo. 
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A stationary Q-lump is a critical point of the energy functional E[tjj], subject to 
Q having some fixed value. The usual (Derrick) scaling argument shows that any 
stationary Q-lump must satisfy 

(2 - n)E q - nE p + nE k = 0. (13) 

For simplicity in what follows, let us take n > 3. Define m > by U"(0) = m 2 ; then, 
near spatial infinity, the Euler-Lagrange equations give \7 2 tp — (m 2 — oj 2 )ij) = 0. So 
in order to satisfy the boundary condition tp — > as r — > oo, we need u < m. 

It is clear from (|13|) that if U > ^m 2 ]^] 2 everywhere, then there can be no 
solution. So K = min [2U(\i/j\)/\i>\ 2 ] has to satisfy K < m 2 . Also, we have 

E p = Ju>^KI= {K/uj 2 )E k > (K/lu 2 )E p , (14) 

where the final inequality comes from ()13|) . As a consequence, we see that uj 2 is 
restricted to the range 

K<uj 2 <m 2 . (15) 

An example which has been studied in some detail is U(f) = f 2 [l + (1 — f 2 ) 2 ]', here 
m 2 = 4 and K = 2, so the range of frequency for Q-balls in this system is y/2 < uj < 2. 
The dynamics of Q-balls in systems such as these turns out to be quite complicated. 
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